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ABSTRACT 

A high quality Keck spectrum of the Ha line in NGC 4395 reveals symmetric exponential wings, 
f v oc e _l '/ CT , with a ~ 500 km s _1 . The wings extend out to XL 2500 km s _1 from the line core, and 
down to a flux density of ^ 10~ 3 of the peak flux density. Numerical and analytic calculations indicate 
that exponential wings are expected for optically thin, isotropic, thermal electron scattering. Such 
scattering produces exponential wings with a ~ l.l(T e (lnr ( r 1 ) ' 45 , where a e is the electron velocity 
dispersion, and r e is the electron scattering optical depth. The Ha wings in NGC 4395 are well fit 
by an electron scattering model with r e = 0.34, and an electron temperature T e = 1.1 x 10 4 K. Such 
conditions are produced in photoionized gas with an ionization parameter U ~ 0.3, as expected in 
the broad line region (BLR). Similar analysis of the [O III] A5007 line yields r e < 0.01, consistent 
with the lower ionization in the narrow line region. If the electron scattering interpretation is correct, 
there should be a tight correlation between r e and the ionizing flux on time scales shorter than the 
BLR dynamical time, or ~ 1 week for NGC 4395. In contrast, the value of a should remain nearly 
constant on these time scales. Such wings may be discernible in other objects with unusually narrow 
Balmer lines, and they can provide a useful direct probe of T e and r e in the BLR. 
Subject headings: galaxies: active — galaxies: Seyfert — quasars: emission lines — galaxies: individual 
NGC 4395 



1. INTRODUCTION 

Early studies of the broad line profiles in Type 1 Active 
Galactic Nuclei (AGN) suggested that the large observed 
line widths ( £ 3000 km s~ x ) are produced by electron 
scattering within the photoionized gas (Weymann 1970; 
Mathis 1970). With the development of photoionization 
models for the Broad Line Region (BLR), it was how- 
ever realized that the expected electron scattering optical 
depth, T e , within the BLR gas is <C 1, which is too low 
to explain the observed line profiles (Davidson & Netzer 
1979). Later studies of the size of the BLR, based on 
reverberation mappings, indicated that the BLR is more 
compact than previously thought (e.g. Rees, Netzer, & 
Ferland 1989), and thus r e is expected to be higher than 
previously thought, though it is not expected to reach 
r e > 1, required to explain the line profiles purely by 
electron scattering. There is also good evidence now that 
the line broadening in the BLR is mostly due to virialized 
bulk motion (e.g. Peterson & Wandel 2000). However, 
electron scattering by a significantly hotter gas outside 
the BLR gas (T ~ 10 6 - 7 AT vs. T ~ 10 4 K ) remains a 
viable option (Shields & McKee 1981; Kallman & Kro- 
lik 1986; Emmering et al. 1992; Bottorff et al. 1997). 
Evidence for electron scattering of continuum photons in 
AGN is provided by wavelength independent polariza- 
tion in some objects (e.g. Antonucci et al. 1993), but 
there is yet no conclusive detection of electron scattering 
effects on the broad line profiles. 

A possible exception may occur in NGC 1068, where 
spatially resolved spectropolarimetry reveals a broad- 
ened Ha line, which may result from electron scattering 
in warm (T ~ lO 5 ^) gas (Miller, Goodrich & Math- 
ews 1991). A spectropolarimetric search for polarized 
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broad lines in low luminosity AGN (Barth et al. 1999) 
revealed a polarized BLR component in the lowest lumi- 
nosity type I AGN, NGC 4395. Based on that detection, 
Barth et al. suggested that higher S/N spectropolarime- 
try may reveal a broadened Ha component, if the po- 
larization is induced by electron scattering. Here we re- 
port on a high quality Keck observation of the Ha line in 
NGC 4395, which indeed reveals extended symmetric ex- 
ponential wings (though their polarization is unknown). 
As shown below, these wings are well fit by optically thin 
electron scattering by electrons at T e = 1.14 x 10 4 AT, with 
T e = 0.34, which are plausible conditions for the BLR 
gas. 

The paper is structured as follows. Section 2 describes 
the observations which revealed the exponential wings in 
the Ha profile of NGC 4395. Section 3 describes the nu- 
merical calculation of the electron scattering line profiles, 
which yield exponential wings, /„ = foe ~ v l a . A fitting 
formula for <t as a function of T e and r e is provided, to- 
gether with a simple analytic derivation. In section 4 
the observed Ha profile is fit with the electron scattering 
model, yielding T e and r e , the origin of the scattering gas 
is discussed, and some predictions are made. The main 
conclusions are summarized in section 5. 

2. OBSERVATIONS 

NGC 4395 was observed with the Keck-II 10 m tele- 
scope using the Echellette Spectrograph and Imager 
(ESI; Sheinis et al. 2002); a slit width of 0".5, and a 
spectral resolution of A/AA = 8,000, or 0.820 A near 
Ha, sampled at 0.259 A per pixel. Two 20 minute expo- 
sures were made on 2002 December 2 UT, yielding a S/N 
ranging from 50 per pixel at the continuum near Ha to 
400 at the line peak (for the reduction details see Laor 
et al. 2006). 
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The continuum level is estimated using a lin- 
ear fit to the observed flux densities of 4.32 x 
10~ 16 erg s _1 cm~ 2 A -1 at observed frame 6480 A, 
and 4.20 x 10~ 16 erg s" 1 cm" 2 A" 1 at 6640 A. The 
implied spectral slope of f\ oc A -115 between 6480 A 
and 6640 A is consistent with the global optical con- 
tinuum spectral slope in the ESI spectrum, although 
the flux at 6640 A appears to be somewhat blended 
with the extended blue wing of the He I A6678.15 line. 
We estimate the continuum placement uncertainty as 
~ 3 x 10~ 18 erg s _1 cm~ 2 A -1 . This uncertainty im- 
plies a potential systematic error of ~ 10% in the net 
line flux at |u| = 2500 km s -1 , increasing to <~ 100% 
at \v\ = 3000 km s _1 . We therefore consider the Ha 
line profile reliable at \v\ < 2500 km s _1 , where any 
systematic flux error is < 10%. The continuum 1 a 
noise level per pixel is 9.5 x 10~ 18 erg s _1 cm -2 , and 
thus the statistical uncertainty per pixel is ~ 30% at 
\v\ = 2500 km s -1 . We find a S/N> 10 per resolution 
element at \v\ < 2000 km s -1 , increasing to > 70 at 
H < 1000 km s _1 . Overall, the high quality of the ESI 
spectrum allows us to reliably probe the Ha profile over 
an unprecedented range of a factor of 1000 in flux den- 
sity, from 3.58 x 10 -14 erg s _1 cm~ 2 A -1 at v — km s _1 
to ~ 3 x 10~ 17 erg s" 1 cm" 2 A" 1 at \v\ = 2500 km s _1 . 

Figure 1 , upper panel, presents the overall Ha line pro- 
file at \v\ < 3000 km s _1 . The nearly pure linear shape in 
log Je versus v is very prominent, implying exponential 
wings. The middle panel demonstrates the asymmetry 
present in the \v\ < 1000 km s _1 region of the broad line 
profile, with excess blue wing flux. Note that the wings 
at \v\ > 1000 km s _1 show an opposite asymmetry (see 
upper panel) , however there is a very good match of the 
red and blue wings when the spectrum is reflected with 
respect to the v = 50 km s _1 point (lower panel). The 
two wings are well matched by an exponential profile, 
Je oc e~ v / a ', with a — 500 km s _1 . These pure symmet- 
ric exponential wings suggest a scattering origin. Below 
we show that exponential wings are indeed a generic pre- 
diction of optically thin isotropic electron scattering. 

3. CALCULATIONS 

3.1. Numerical 

A complete calculation of electron scattering line pro- 
files involves a calculation of the photon diffusion simul- 
taneously in both real space and in energy space. Such 
a calculation can be carried out using a Monte Carlo ap- 
proach, which is straightforward to implement, though 
the finite number of photons used tends to produce non 
smooth curves (e.g. Lee 1999) compared to analytic solu- 
tions. Analytic solutions for the local diffusion in space, 
described by the radiative transfer equation, and diffu- 
sion in energy, described by the energy redistribution 
function, have been derived and solved numerically us- 
ing various simplifying assumptions (e.g. Chandrasekhar 
1960; Mihalas 1978; and references therein). Here we 
make two major simplifying assumptions. First, the 
electron scattering gas is assumed to be optically thin, 
thus a detailed radiative transfer solution is not required. 
This simplification is justified by the observation that the 
dominant effect on the line broadening in the BLR ap- 
pears to be gravity dominated bulk motion, rather than 
electron scattering. Gravity is indicated by the good cor- 
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Fig. 1. — Upper panel: The observed Ha profile in NGC 4395 
(solid histogram), and the reflected profile (dotted histogram). 
Note that the blue wing is higher than the red wing at \v\ < 
1000 krns" 1 , while the reverse is true at \v\ > 1000 km s 1 . Middle 
panel: An expanded view of the line core. The narrow Ha line dom- 
inates at \v\ < 150 km s _1 , and the narrow [N II] lines dominate 
around v ~ 950 km s _1 , and v ~ —670 km s _1 . Lower panel: An 
expanded view of the extended wings at 1000 < \v\ < 3000 km s _1 . 
The velocity scale in this panel was shifted from v to v — 50 km s _1 . 
Note the remarkable symmetry of the wings in the shifted scale, 
and their nearly pure exponential form. The solid line is an expo- 
nential function with a = 500 km s _1 . 



respondence between black hole mass estimates assum- 
ing a virializcd BLR and mass estimates based on the 
host bulge properties (Laor 1998; Ferrarese et al. 2001). 
Second, we assume an isotropic electron velocity distri- 
bution, and an isotropic illumination of the scattering 
electrons, which is likely to be a good approximation for 
the integrated scattered radiation. These assumptions 
lead to the following remarkably simple energy redistri- 
bution function, derived by Rybicki & Lightman (1979, 
hereafter RL79. Eq. 7.24 there), 

j(x,P) = ^^[(l + (3)x-(l-P)} fori T |< a; <l 

(1) 

j{^P)= l -^[(l + P)-x{l-m forl< a; <i±| ) 

(2) 

j(x,/3)=0 for x > YZTJ) or x < Y+~^> ( 3 ) 

where x = e/eo, and eo,e are the observed photon en- 
ergy before and after scattering, v = [3c is the scattering 
electron velocity, and j(x,(3) is the number of photons 
scattered per dx, which obeys the photon conservation 
Jq° j{ x > P)dx = 1. The expressions above for j{x,(3) ap- 
ply for elastic scattering in the electron rest frame, which 
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Fig. 2. — Upper panel: The single velocity redistribution function 
j(x,/3) for T e = 10 4 K (triangle shaped curve), and the thermal 
averaged function jt(x), scaled to match j(x,/3) at x = 1 (i.e. at 
v = km s _1 ). Note that despite the thermal averaging the core 
of the redistribution function is almost unchanged. Lower panel: 
The effects of repeated scattering. The narrow Gaussian is the 
unscattered line profile, and the broader curves correspond to /" 
with n = 1 — 10. The number of photons (area under each curve) 
is assumed to be conserved, and the unscattered line profile was 
scaled to match the peak of the n = 1 curve. 



is a very good approximation for the applications below, 
where -C 1 and hu <C m e c 2 . 

To obtain the redistribution function for a thermal 
electron distribution, jt(x), one needs to integrate the 
above expression for j(x,0) over a thermal velocity dis- 
tribution function. Thus 

A /*°° 

j t (x) = —— j(x,(3)v 2 e-WV 2 dv (4) 
b V 7r Jo 

where b = 2kT e /m e , and T e is the electron temperature. 
Figure 2 (upper panel) shows j(x,0) and jt(x). Note 
that despite the thermal averaging, j t (x) remains sharply 
peaked. 

Since the scattering medium is optically thin, each 
scattering redistributes the observed line profile through 
jt(x), and reduces the number of photons by a factor 
of T e . Thus, the line profile following n + 1 scattering is 
given by a convolution of jt(x) with the n times scattered 
profile, i.e. 

fl n+1) =r e f:,j t [(v-v')/c]dv' (5) 
Jo 

where /" is the profile of a line scattered n times. Since 
<C 1, a simple Doppler relation de/e — dv/c can be 
assumed. Figure 2 (lower panel) shows /°, an arbitrary 
seed line profile, and /" for n = 1 — 10, illustrating how 
the line gets broader and shallower with repeated scatter- 
ing by a given electron population. The total scattered 




-4000 -2000 2000 

v (km s -1 ) 



Fig. 3. — The dependence of /^ cat on T e and r e . Upper panel: 
The dependence on T e for r e = 0.3. The line wings approach a 
linear shape in log /|j cat versus v, i.e. an exponential form /° cat oc 
e~ v >' 7 , at high v. Note that cr increases significantly with T e . 
Lower panel: The dependence on r e for T e = 2 X 10 4 K. The 
integrated scattered flux is proportional to r e , however a is only 
weakly dependent on r e . 



line profile is then given by 

oo 

71=1 

Figure 3 shows the results of numerical calculations of 
jscat £ Qr a ran g e G f y e anc j Te va i ues . As both T e and 

r e increase, the line profile gets broader, with a rela- 
tively weak dependence of the width on r e and a stronger 
dependence on T e . Note that in all cases log/^ cat ap- 
proaches a linear function of v in the extended wings, or 
cquivalently 

/r' - he- v/ °. (7) 

Figure 4 shows a plot of a, measured from the numerical 
calculations, as a function of r e and T e . The solid lines 
is a fit to the numerical results. The fitting function has 
the simple form, 

( 7 = l.l<7 e (lnr- 1 )- - 45 , (8) 
where a e is the thermal electron velocity dispersion, 

cr e = ^JkT e /m e = 389.2T 4 1/2 km s~\ (9) 

and T 4 = T e /10 4 . This fitting function is accurate to 
1.5% in the range 0.01 < r e < 0.8. 

3.2. Analytic Derivations 

3.2.1. An Oversimplified Derivation 

One may attempt to understand the above fit equa- 
tion based on the following very simple derivation. This 
derivation is actually oversimplified, and the resulting 
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expression for a lacks accuracy, but it provides insight 
for understanding the basic reason for the formation of 
exponential profile wings. 

This derivation follows closely the derivation of the 
power law spectral form produced through optically thin 
electron scattering by relativistic electrons, as described 
in RL79 (section 7.5 there). Here we have optically thin 
electron scattering by non relativistic electrons. Let a e 
be the mean electron velocity. The energy amplification 
of a photon scattered at right angle is 

ei/e ~ l + cr e /c. (10) 

After n such scattering the photon energy will be 

e„ ~ e (l + cr e /c)™, or e„ ~ e (l + na e /c), (11) 

where the last step is valid for na e /c <C 1. Changing to 
velocity scale 

v n = (e„/e - l)c, gives v n ~ na e . (12) 

The fraction of photons scattered n times is r™, and thus 



/K)~/(0)r™ or f(v n )~f(0)r^ 
This is equivalent to 

/W~/(0)e b W", 
i.e. an exponential line profile, where 

/(«) oc /(0)e" , ' /<7 and <j = a e /\nT~ 



(13) 



(14) 



(15) 

The crucial difference here from relativistic electron scat- 
tering is that the energy increment with repeated scat- 
tering is a constant linear factor, rather than a constant 
multiplicative factor, which would lead to a power law 
rather than an exponential form. 

The above derivation of a does not yield the correct 
dependence on In r e . The oversimplification in the above 
derivation is in the assumption that the energy redis- 
tribution function can be approximated as a 6 function 
at ei/eo — 1 + <J e /c, while as Eqs.1-3 above show, the 
photons actually get redistributed nearly symmetrically 
from 1 — a e /c to l + a e /c. For highly relativistic electrons 
the strong Doppler beaming produces a highly asym- 
metric and strongly peaked redistribution function (at 
ei/eo — 7 2 ), and the simplified derivation of the power 
law spectral slope through repeated scattering (RL79) is 
adequate. 

3.2.2. A Simplified Derivation 

Taking into account the shape of the redistribution 
function leads to a more accurate dependence of a on 
r e , as described below. 

The thermal redistribution function can be approxi- 
mated as a Gaussian with the velocity dispersion of the 
scattering electrons, jt(x) oc e~ v l h . The seed line pro- 
file is assumed to have a smaller velocity width than that 
of the scattering electrons, so it can be conveniently ap- 
proximated as a S function, = 5(v). Thus, the first 
scattering (Eq.5) yields = jt(x), and repeated scat- 
tering amount to repeated convolutions with jt(x). Such 
convolutions yield Gaussians with a velocity dispersion 
which grows as y/n, i.e. we can use the approximation 
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(16) 




Fig. 4. — The dependence of a on r e for two values of T e . 
The filled triangles are values measured from numerical calcula- 
tions of /5 cat , and the solid lines are analytic fits of the form 
a = l.l<r e (ln r^ 1 ) -0,45 to the numerical results (see text). This fit 
is accurate to 1.5% in the displayed range. 



where K is an arbitrary constant, and the 1/t/ti term 
comes from the Gaussian normalization. The total scat- 
tered line profile is then given by 



rscat jv- \ ^ e 



-v 2 /nb 2 



(17) 
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Unfortunately, there is no simple analytic solution for the 
sum of this series. Note, however, that only a few terms 
in the above series contribute at any given v. The contri- 
bution of terms with n <C v 2 /b 2 is small because of the 
large negative exponential factor, while the contribution 
of terms with n 3> v 2 /b 2 is reduced by the r™ / \/n factor. 
The maximum contribution to /^ cat in the above sum is 
provided by the term n max , which is obtained from the 
requirement 

— — = (hiT e 

an In 

which holds at n 



(18) 



quadratic equation for n ma 



n 2 b 2 ' 

Since /" ^ one gets a 



2 In r e b 2 In r e 



= 0. 



(19) 



The second term can be neglected compared to the first 
and third terms when n > 1 and r e < 0.6. Thus, the 
solution is 

v 

Wmax = , (20) 

6vhiT e 

in the general expression for /" gives 



Substituting n n 



^v/ibV^Te 1 ) 



(21) 
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We further show that /" indeed has a narrow peak 
around n — n max . We define the peak width, An, 
through the relation /" max+An //" max _ q wn j c h yields 
after some steps the condition 

r 2(A„) 2 /" ^ 5) (22) 

or 

An ~ Otin^/VlnTe 1 , (23) 

i.e. only the nearest 1-2 bins to n max contribute signifi- 
cantly for n ^ 10. Thus, at each v one can approximate 
the above sum over n by just a single term at n = n max , 
i.e. /^ cat = /™ max . To verify that /^ cat has the correct 
exponential form we calculate 

1 = rfln/f at = 2y/j^g _ 1_ 

Thus, a pure exponential line profile, i.e. a independent 
of v, is obtained when the second term on the right is 

negligible, which occurs at v > The expo- 

nential slope is a = 6/2ylnTe , or a = 0.71a e 
which matches well the lnr e dependence of the fitting 
function to the numerical solution (Eq.8). 

4. RESULTS AND DISCUSSION 
4.1. The best jit electron scattering model 

To model the shape of the electron scattering wings 
one needs an estimate of the seed line profile. As shown 
below, there is no evidence for detectable electron scat- 
tering wings for lines originating in the narrow line region 
(NLR), and we therefore use only the broad component 
of Ha for the seed profile. To obtain the BLR Ha compo- 
nent, we subtract off the narrow Ha component and the 
narrow [N II] lines from the observed Ha profile. The 
broad component underlying the narrow lines is inter- 
polated using a spline fit, where we use the minimum 
number of spline points required to get both an accept- 
able x 2 j and also matching profiles for the narrow [N II] 
doublet lines (Laor et al. 2006). 

The observed broad line core at \v\ < 1000 km s _1 
shows a noticeable asymmetry, with up to ~ 50% higher 
flux in the blue wing (Fig.l, middle panel). This ap- 
pears to be a persistent effect in NGC 4395, as it is also 
present in high resolution Keck spectra taken about a 
year and eight years earlier (Fig.l in Laor et al. 2006). 
In contrast, the wings at \v\ > 1000 km s _1 are generally 
symmetric to < 10% (with respect to v — 50 km s _1 , 
Fig.l, lower panel). Since our electron scattering model 
generates symmetric line profiles, one must assume that 
the observed asymmetric core profile represents the un- 
scattered seed line. To estimate the shape of the seed 
line wings, the observed profile is multiplied by a weight 
function of the arbitrary functional form e~ v l w , which 
produces a smooth cutoff in the the wings at \v\ > w. 
We find that w = 600 km s _1 at v > km s _1 , and 
w = 800 km s _1 at v < km s _1 , are about the min- 
imum values required to exclude the inner asymmetric 
core emission, and allow the electron scattering model to 
fit only the extended symmetric wings. The exact value 
of w cannot be accurately constrained, and the true seed 
profile may be broader than assumed here. The seed pro- 
file is further multiplied by e~ r<! , which represents the 



wavelength independent suppression due to the electron 
scattering, where t c is determined below. 

The seed line profile obtained above is now used as 
the input line profile for our electron scattering code. 
The model free parameters are T e and r e , and we iterate 
over different values for these parameters until a satis- 
factory fit is obtained between the observed broad Ha 
profile, and the sum of the scattered + seed line pro- 
files (see Figure 5). The scattering medium is assumed 
to have a bulk velocity of v — 50 km s _1 , with respect 
to the narrow line peak, as suggested by the observed 
asymmetry of the extended wings. The best fit is ob- 
tained for T e = 1.14 x 10 4 K, and r e = 0.34, which yields 
X 2 = 282 for 243 degrees of freedom, where \ 2 1S cal- 
culated for \v\ < 2500 km s _1 , excluding the line core 
region, —900 < v < 1200 km s _1 , affected by the nar- 
row lines. The best fit T e and T e depend on the assumed 
shape for the seed line profile. Increasing the wings cut- 
off velocity, w, reduces the best fit r e and increases T e . 
For example, adopting w = 1200 km s _1 yields a best fit 
solution with T e = 1.7 x 10 4 K, and r e = 0.17. However, 
at such a high w the seed line contribution becomes neg- 
ligible (< 20% of total line flux) at \v\ > 1500 km s" 1 , 
compared to \v\ ^ 1000 km s _1 for the minimum w so- 
lution. The observed wings are well fit by a featureless 
exponential profile already at \v\ ^ 1000 km s _1 (Fig.l, 
lower panel), which suggests this region is dominated by 
a single component, as we find for the minimum w solu- 
tion. We therefore adopt below the best fit minimum w 
solution, although one cannot rule out a higher w solu- 
tion, where the seed line profile and the scattered wings 
conspire to produce a single featureless exponential pro- 
file at \v\ > 1000 km s" 1 . 

Once the seed line profile is set, there is essentially no 
degeneracy in the solution, as r e is fixed uniquely by the 
ratio of seed to scattered photon flux ratio, and T e is 
then fixed by r e and the observed exponential slope, a. 
The formal uncertainties in T e and r e due to the avail- 
able S/N here are < 3 — 5%. However, our assumption 
of an isothermal and isotropic scattering medium may 
be a significant oversimplification, and systematic errors 
£ 10% are quite possible due to our simplified modeling. 

4.2. Where does the scattering takes place? 

The best fit T e is typical of photoionized gas in the 
BLR or the NLR. To clarify whether the scattering gas 
resides in the BLR or in the NLR, we looked for evi- 
dence for electron scattering wings in the profile of the 
[O III] A5007 line, the strongest narrow line in the spec- 
trum of NGC 4395. Figure 6 shows the continuum sub- 
tracted line profile. Weak broad wings appear to be 
present at |u| > 500 km s _1 . To determine the implied r e , 
one needs an estimate of the seed line profile. As for the 
Ha line, we multiplied the observed profile by a weight 
function of the form e~ v l h , using here b — 400 km s _1 . 
We fixed T e at 1.14 x 10 4 K, the value measured for Ha, 
and obtained a reasonable fit to the wing profile using 
r e = 0.01. However, weak broad wings are also present 
in the instrumental line spread function (LSF) at a level 
of ~ 10~ 3 of the peak line flux density. Convolving the 
seed [O III] A5007 line profile with the instrumental LSF 
produces wings which are also roughly consistent with 
the observed wing flux (Fig. 6). We thus conclude that 
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Fig. 5. — The observed Ha profile (solid histogram), and the 
best fit optically thin electron scattering model (solid line). The 
input seed line profile and the electron scattered profile are also 
shown. The best fit parameters and the fit \ 2 are indicated on the 
plot. 




Fig. 6.— The observed [O III] A5007 line profile (solid his- 
togram), the assumed seed line profile (dotted histogram)), and 
the best fit optically thin electron scattering model (solid line) 
having r e = 0.01 and T e fixed at 1.14 X 10 4 K. The wings are also 
consistent with instrumental broadening (dashed line), and there- 
fore we deduce that r e < 0.01. The feature at v ~ 500 km s _1 is 
the He I A5015.68 line. 



most, or possibly all, of the [O III] A5007 wing flux may 
be instrumental, and therefore one can only put an upper 
limit of r e < 0.01 on a T e ~ 10 4 K scattering medium at, 
or outside, the NLR. 

In contrast to the [O III] A5007 line, the instrumental 
LSF has a negligible effect on the Ha line wings. For 
example, the instrumental LSF, acting on the narrow Ha 
component, produces a flux density at |u| = 1500 km s _1 
which is < 10~ 4 of the peak flux density (Fig. 6), while 
the observed Ha flux density at \v\ — 1500 km s _1 is 
~ 5 x 10~ 3 of the peak flux density (Fig. 5). In addition, 
the instrumental LSF scatters only ~ 1% of the line flux 
into broad wings, while the observed broad wings flux 



in Ha amounts to 34% of the line flux, again indicating 
that the exponential Ha wings cannot be an instrumental 
effect. 

Since the [O III] A5007 line shows no scattering wings, 
the scattering wings of Ha must originate in gas at radii 
smaller than the NLR, most plausibly within the pho- 
toionized gas which produces the BLR emission lines. 
An optical depth r e = 0.34 implies an electron column 
of 5.1 x 10 23 cm~ 2 , or cquivalcntly a minimum total H 
column of about 4.6 x 10 23 cm~ 2 . Photoionized gas typ- 
ically has an ionized, H II, surface layer with a column 
of S; on :=y 10 23 U cm~ 2 , where U = n 7 /n e is the ioniza- 
tion parameter, and n 7 and n e are (respectively) the H 
ionizing photon and the electron number densities. To 
yield the required r e within this surface layer requires 
U ~ 5. This value of U is much larger than the typical 
U ~ 0.1 at the BLR, and it implies T e ~ 5 x 10 4 K within 
the H II layer. Such a high T e is clearly excluded by our 
measurement of T e = 1.14 x 10 4 K within the Ha scatter- 
ing layer. However, below the surface H II layer lies an 
extended layer of partially ionized gas, where the lower 
ionization lines (Mg II, Fe II), and most of the Balmcr 
lines originate. The Ha scattering wings may therefore 
originate in this deeper layer. 

The ionization structure of the partially ionized layer 
cannot be deduced by simple arguments, as made above 
for the H II layer. We therefore use the photoionization 
code CLOUDY (Fcrland et al. 1998) to calculate r e , as 
a function of U and total column Eh, for models with 
a significant H I column. The partially ionized layer is 
heated by the harder, more penetrating radiation. We 
assume a spectral energy distribution (SED) having an 
energy spectral slope of — 1 from 1 eV to 10 eV, —1.5 from 
10 eV to 1 keV, -0.2 from 1 keV to 40 keV, and a cutoff 
above 40 keV, which is consistent with the SED mea- 
sured for NGC 4395 by Moran et al. (1999). This X-ray 
hard SED produces about twice the r e obtained for the 
softer Mathewes & Ferland (1987) SED, commonly used 
for photoionization calculations. The observed r e = 0.34 
is obtained for S H = 10 24 cm~ 2 and U — 0.3, well within 
the plausible range for the BLR. The CLOUDY model 
calculations show there is a significant temperature gra- 
dient within the electron scattering layer. The surface 
temperature is 3.4 x 10 4 K, dropping to 1.6 x 10 4 K at 
a depth where r e = 0.07, 7.2 x 10 3 K at r e = 0.19, and 
7 x 10 3 K at r e = 0.3. Calculating the electron scattering 
profile of a nonisothermal slab requires a complete radia- 
tion transfer calculation, including the Ha emissivity as 
a function of depth, which is well beyond the scope of 
this paper. However, our isothermal electron scattering 
model may capture the r e weighted mean temperature of 
the electron scattering medium, which is rather close to 
the T e ~ 1.1 x 10 4 K derived above. 

4.3. Is it really electron scattering? 

Although our best fit model parameters, T e and r e , are 
tightly constrained by the observed Ha profile, it is not 
clear that this model is the only possible interpretation. 
One may assume that no electron scattering is present, 
and that the Ha gas emissivity as a function of velocity 
in the BLR just happens to produce an exponential pro- 
file. How can the electron scattering model be tested? A 
clear prediction of our model is that r e should be tightly 
correlated with the ionizing flux. As the ionizing flux 
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increases the column of ionized gas must increase, and 
the relative fraction of line flux in the scattering wings 
should increase. If the BLR gas has a large column, so it 
is ionization bounded rather than matter bounded, then 
the relation between r e and the ionizing flux should be 
linear. In contrast, the slope of the exponential wings 
a should remain nearly constant. This occurs since T e 
is only weakly dependent on U (for U < 1), and thus 
should vary only slightly. In addition, although r e is 
expected to vary linearly with U, a is only weakly de- 
pendent on r e , and would thus not be significantly af- 
fected. The ionization structure responds on a timescale 
of ~ 4 x 10 12 /n e seconds, or ~ 400 seconds for the BLR, 
which is much smaller than the BLR light crossing time 
of ~ 1 hour (Peterson et al. 2005). Thus, one expects the 
measured r e and T e to respond to changes in the ionizing 
continuum on timcscales as short as ~ 1 hour. However, 
since the dynamical time in the BLR of NGC 4395 is 
<~ 1 week, structural changes may occur on this time 
scale, or longer, which could eliminate the correlation of 
T e and T e with the ionizing flux (e.g. due to changes in 
the BLR density and column density). 

Indirect evidence for electron scattering in NGC 4395 
is provided by the spectropolarimetry of Barth et al. 
(1999), who find a wavelength independent polarization 
of 0.67 ± 0.03%. The similar low polarization of the 
continuum and the broad lines suggests the scattering 
most likely occurs by a low T e scattering medium outside 
the BLR, rather than within the BLR gas, as suggested 
here 1 . Barth et al. suggested that "A clear sign of 
electron scattering would be a broadened Ha profile in 
Stokes flux," . The exponential Ha wings revealed in this 
study are probably not related to the Barth et al. pre- 
diction, unless scattering off the BLR clouds can produce 
the observed continuum and line polarization (Korista & 
Ferland 1998). Spectropolarimetry at a higher S/N and 
higher spectral resolution is required to better constrain 
the location and identity of the polarizing medium in 
NGC 4395. 

Indirect support for the electron-scattering-within-the- 
BLR interpretation is provided by the fact that best 
fit values of the two free model parameters, T e and 
T e , can apparently be obtained using a single parame- 
ter, U ~ 0.3. This value of U is consistent with the 
range of U values typically deduced in AGN, though 
it is significantly higher than deduced by Kraemer et 
al. (1999) for this object. Can the required U and Eh 
be accommodated within the small BLR of NGC 4395? 
To estimate the physical size of such a BLR compo- 
nent we use the mean "visit 3" flux density of fx = 
3 x 10~ 15 erg" 1 s" 1 cm~ 2 A" 1 at 1350A, taken from 
Table 3 in Peterson et al. (2005, note an Erratum with 
corrected flux), which matches well an extrapolation of 
our ESI spectrum to the UV. At a distance of 4.3 Mpc 
(Thim et al. 2004), the corresponding monochromatic 
luminosity is \L\ = 1.05 x 10 40 erg" 1 s" 1 , corrected for 
Galactic extinction of 0.175 mag at 1350A (Schlegel et 
al. 1998 and Table 1 in Mathis 1990). Assuming the 
ionizing luminosity is ~ 3AL> at 1350A, and a mean ion- 
izing photon energy of ~ 2 Rydberg, the implied ionizing 



photon flux from the source is N ~ 7 x 10 50 s _1 . Us- 
ing U = N I (inr 2 cn e ) we get U ~ 20r 14 2 n^ 1 , where the 
distance of the emitting gas is r — 10 14 ri4 cm, and its 
density is n e — 10 10 ni cm~ 3 . The measured value of 
U = 0.3 implies a distance 7*14 = 8n^ ' 5 . The thickness 
of the gas slab is du = £24/^10 cm j an d thus for £24 = 1 
found here we get d/r = 0.12n^ ' 5 . Therefore, d/r 0.1 
for n e ^ 10 10 cm~ 3 , and the implied column of Ha emit- 
ting gas can be accommodated within the compact BLR 
of NGC 4395, for reasonable values of n e . 

Why then are such exponential wings not commonly 
seen in other AGN? NGC 4395 is unique among type 
I AGN in having extremely narrow Balmer lines. The 
broad Ha FWHM here is just 520 km s _1 , while the 
typical Ha FWHM in AGN is ~ 3000 km s _1 , with a 
few times broader line base. Thus, the electron scattering 
effect of the BLR gas will have a negligible effect on the 
emission line width in most AGN. However, one can make 
a clear prediction that the exponential wings should be 
discernable in other objects with unusually narrow lines. 
A reliable detection of such wings requires a high S/N 
spectrum (e.g. 10-70 here at 1000 < |u| < 2500 km s" 1 ). 
For example, the H/3 line wings of NGC 4395, and the 
Ha line wings in a Keck spectrum of POX 52 (Barth 
et al. 2004, kindly provided by A. Barth), are broadly 
consistent with an exponential form, but the available 
S/N in both cases is too low to constrain the profile shape 
as accurately as done here. We note in passing that as 
long as r e remains low ( £ 0.5), electron scattering would 
not have a significant effect on measurements of the line 
FWHM, and thus on estimates of the black hole mass 
(e.g. Kaspi et al. 2005). 

5. CONCLUSIONS 

High quality Keck observations of the Ha line in 
NGC 4395 reveal symmetric exponential wings of the 
form f v oc e~ v l a , with a ~ 500 km s" 1 . The wings are 
fit well by an isothermal and isotropic electron scattering 
model with r e = 0.34, and T e = 1.14 x 10 4 K, plausibly 
generated in photoionized gas with U ~ 0.3. The lack of 
electron scattering wings for the [O III] A5007 line, and 
the value of U indicate that the scattering occurs within 
the BLR gas where Ha is formed. 

The electron scattering interpretation can be tested by 
looking for a strong correlation between r e and the ion- 
izing flux. High quality spectra of other AGN with very 
narrow Balmer lines will provide a further test of the 
electron scattering interpretation. If the electron scat- 
tering origin is verified, then the extended Ha wings can 
provide us with a new direct probe of T e and T e within 
the BLR gas, and their time variability. 
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1 Inspection of the 2-D spectral images of NGC 4395 obtained nuclear light 

with the HST STIS suggests some spatially extended scattered 
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